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Abstract 



We analyze the solar and the atmospheric neutrino problems in the context of 
three flavour neutrino oscillations. We assume a mass hierarchy in the vacuum 

(N ■ 

' mass eigenvalues ^ ^> /i| > A 1 ?, but make no approximation regarding the 

in 

<3\ ' magnitudes of the mixing angles. We find that there are small but continuous 

Q_i' bands in the parameter space where the constraints imposed by the current 

measurements of 71 Ga, 37 Cl and Kamiokande experiments are satisfied at la 
level. The allowed parameter space increases dramatically if the error bars 
^ ' are enlarged to 1.6a. The electron neutrino survival probability has different 



energy dependence in different regions of the parameter space. Measurement 
of the recoil electron energy spectrum in detectors that use v— e scattering may 
distinguish between some of the allowed regions of parameter space. Finally 
we use the results for the parameter space admitted by the solar neutrinos as 
an input for the atmospheric neutrino problem and show that there exists a 
substantial region of parameter space in which both problems can be solved. 

PACS numbers: 14.60. Gh, 96.60.Kx, 95.30.Cq, 96.40.Tv 
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I. INTRODUCTION 



The solar neutrino problem has been an interesting and intriguing phenomenon in neu- 
trino physics for a long time. The different solar neutrino experiments observe differing 
fractions of the neutrino flux predicted by the Standard Solar Model (SSM) [fl],EJ. The 
oldest of the solar neutrino experiments is the 37 C7 experiment at Homestake. Its energy 
threshold is 0.814 MeV and it can detect the neutrinos from 7 Be (E u = 0.862 MeV) and 8 B 
[E v < 14.02 MeV) reactions. In the standard solar model (SSM) of Bahcall-Pinnsonneault 
0, the capture rate in the 37 Cl experiment is predicted to be 8.0 ± 1.0 SNU. However, the 
measured rate is only [|J 

R C i = 2.55 ± 0.25 SNU. (1) 

The water Cerenkov detector at Kamioka, with a threshold of 7.5 MeV, can detect only the 
neutrinos from the upper end of 8 -B spectrum and the Kamioka result f| is, 

y Kam = * Kam =0-51 ±0.07, (2) 

J^Kam-.SSM 

which is the ratio of the observed neutrino flux to that predicted by the SSM. The gallium 
experiments SAGE and GALLEX, with energy threshold of 0.233 MeV can detect the neu- 
trinos coming from the dominant p — p reaction (E v < 0.42 MeV) as well as the neutrinos 
from 7 Be and 8 B reactions. Their measured rates are |||7| 

Rsage = 69 ±11 ±6 SNU 
Rgallex = 79 ± 10 ± 6 SNU 



and the average is 



R G a;avg = 74 ± 8 SNU (3) 



as opposed to the SSM prediction of 131.5 SNU. 

A rough model independent analysis of these results indicates that the low energy neutri- 
nos from the p — p reaction suffer very little suppression whereas the higher energy neutrinos 



are suppressed to a large extent Recently it was pointed out that if neutrinos have no 

properties beyond those in the standard electro- weak model (i.e. if they are massless), the 
measurement of Kamiokande, together with that of 37 CI experiment, implies that the 7 Be 
neutrinos must be suppressed by more than 90% [0,0] ■ 

In addition there exists an anomaly in the ratio of observed muon neutrinos to electron 
neutrinos in the earth's atmosphere. These neutrinos are produced from the decay of 7r ± 
and which are in turn produced by cosmic rays interacting with the atmosphere. The 
ratio is roughly two as suggested by the Monte-Carlo calculations whereas both Kamioka 



11] and 1MB |L2| report that the ratio is only about half of that predicted by the Monte- 



Carlo calculations. The results for this ratio are also available from three other groups 



using the tracking detectors, namely the NUSEX |R|, Frejus |14| and SOUDAN-II fl5 
collaborations. The data from the NUSEX collaboration seems to be in agreement with no- 
anomaly situation. Similar conclusion is obtained from Frejus data if all the contained events 
are considered. However, if only fully contained results are taken into consideration, there is 
a suppression. The SOUDAN-II results are consistent with the results obtained with water 
Cerenkov detectors. It should be noted that the statistics in the tracking experiments is not 
as high as the water Cerenkov experiments. Evidently any solution of the solar neutrino 
puzzle must incorporate simultaneously a solution of the atmospheric neutrino problem Jl6| . 



A satisfactory solution to the solar neutrino problem should be able to explain not only 
the total deficit that is observed but the differential suppression observed at low and high 
energies. Solutions based on astrophysics or nuclear physics ascribe the deficit to smaller 
solar core temparature or smaller cross sections for the nuclear reactions talking place in 
the sun. Recent model independent analyses suggest that these solutions cannot describe 
the results of 37 CI and Kamiokande simultaneously Particle physics based solutions 

attempt to account for the deficit by assuming that the neutrinos have interactions beyond 
those of the standard electro-weak model. If the neutrinos possess small mass, an electron 
neutrino can oscillate into a neutrino of another flavour [|17||. The amplitude of oscillation 



is a function of the mass squared differences, the mixing angles between neutrino flavours 



and the neutrino energy. If one of the mass square differences is of the order of the effective 
mass squared arising from v e — e interaction, the matter effects can enhance the mixing to 
its maximal value and the amplitude for u e oscillating into another flavour will be very large 
0. This is the so called MSW effect. 

Matter-enhanced oscillations have been studied thoroughly in the scenario where only two 
flavours, u e and u^, mix with each other |l9^pT|. The vacuum oscillation here is controlled 



by the two parameters, the mass square difference S21 = m\ — m\ and the mixing angle u. 
Matter effect is taken into account by adding to the mass squared of u e , the term 

A{r) = V2G F n e (r) x 2E, (4) 

which is proportional to the electron number density in the Sun n e (r), where r is the radial 
distance from the centre of the Sun. The maximum value of A occurs at the core and is 
roughly 10~ 5 E eV 2 , where E is the neutrino energy in MeV. The mixing angle u m in the 
presence of matter is given by, 

5 2 icos2u-A 

cos2u m = . (b) 

'(5 2 i cos 2u - A) 2 + (5 21 sin 2c;) 2 



The MSW resonance condition is, 

A = 621 cos2cj. (6) 

Note that, if the resonance condition is to be satisfied, A core > 5 2 iCos2co>, which implies 
that u m > 7r/4 at core. At resonance it becomes 7r/4 and approaches its vacuum value after 
passing through the resonance. 

The probability for an electron neutrino produced in the solar core to be detected as an 
electron neutrino on earth, averaged over the time of emission and the time of absorption, 
is given by 

(Pee) = COS 2 U COS 2 UJ m + SU1 2 UJ sin 2 UJ m — X\2 cos 2uj cos 2uj m , (7) 

where uo m is to be evaluated at the point of production and x\2 is the probability of a non- 
adiabatic jump between the matter dependent mass eigenstates. If the variation of the solar 
density in the resonance region is slow enough, the adiabatic condition 



is satisfied and the matter dependent mass eignestates evolve adiabatically and there are no 
transitions between them. If (§) is not satisfied, then there will be non-adiabatic transitions 
between the two matter dependent mass eigenstates in the resonance region and the proba- 
bility of this jump has the general form exp(—C/E). The term C has dimensions of energy 
and is some function of 821, to and the derivative of the solar density. The expressions for 
C for various density profiles are tabulated in Ref. |[22| . For linear density variation in the 
resonance region, the jump probability is given by the Landau-Zener formula 



X12 = exp 



7T 

"2 7 . 



(9) 



The predictions for the rates of various experiments are obtained by convoluting the SSM 
neutrino fluxes with the expression for survival probability in (0). A fit to the data from 
71 Ga, 37 Cl and Kamiokande experiments yields solutions in two regions in the #21 — sin 2 2a; 
plane, one with small vacuum mixing and one with large vacuum mixing: 

6 21 ~ 6.1 x l(T 6 eV 2 and sin 2 2co ~ 0.0065, 

S 21 ~ 9.4 x 10" 6 eV 2 and sin 2 2u ~ 0.62. (10) 

In case of the small mixing angle solution, the resonance occurs for neutrinos with energy 
greater than 0.6 MeV. Therefore, the p — p neutrinos (whose maximum energy is 0.42 MeV) 
are unaffected whereas the neutrinos with energy greater than 0.6 MeV are almost com- 
pletely converted into v^. But the measurement of Kamiokande shows that the neutrinos 
with energy greater than 7.5 MeV are suppressed by only a factor of 0.5. This can be ac- 
commodated through the non-adiabatic jump x\ 2 in (|7|). If C ~ lOMeV, or equivalently 
821 sin 2 2uj ~ 4 x 10 -8 eV 2 , then x\ 2 is negligible for energies less than 5 MeV, but becomes 
appreciable at higher energies and (Pee) satisfies Kamiokande constraint. The energy depen- 
dence of (P ee ) in this case is precisely of the form that is required to satisfy the data from the 
three solar neutrino experiments. In the case of the large angle solution, the non-adiabatic 



effects are totally negligible and the (P ee ) is about 0.55 below 0.5 MeV and slowly falls 
to about 0.35 around 5 MeV after which it it remains almost independent of the neutrino 
energy. 

In the case of two flavour oscillations, the area of the parameter space, which can satisfy 
all the three constraints at la level, is very small. Especially, in the case of the small angle 
solution, the requirement that the resonance should occur around 0.6 MeV uniquely fixes 
the value of 621- The requirement that the 7 Be neutrinos should be completely suppressed 
and that the high energy S B neutrinos should have a suppression of about 0.5 determines the 
product S21 sin 2 2oj almost exactly. Therefore, there is very little leeway in the allowed values 
of 821 and sin 2 2uj. An appreciable region of parameter space is allowed only at 95% C.L. ( or 
2.4 a level). In addition, this simple picture is inadequate to simultaneously explain the solar 
and atmospheric neutrino deficits since the mass squared differences required are in vastly 
different regimes. To explain the atmospheric neutrino anomaly on the basis of two-flavour 
vacuum oscillations, one requires a mass squared difference of the order of 10 _1 — 10~ 3 eV 2 , 
with a large mixing angle. This must be compared with the best fit to the data in the case 
of solar neutrino problem given in eq.(10). Therefore one has to necessarily consider the 
scenerio in which all the three neutrinos participate. This of course is also a more realistic 
situation since the LEP experiments have already pinned down the number of light neutrino 
generations to be three. 

Three flavour oscillations were considered previously p3|-p5]|. However, the uncertainties 



in the Gallium experiments have come down significantly in recent times and the parameter 
region allowed by the current data will be much smaller. Recently Joshipura and Krastev 
26| have attempted a complete solution of the solar and atmospheric neutrino problems in 



the three generation frame work. They present a combined analysis of these two problems 
in the framework of the MSW effect and indeed show that there exists a parameter space in 
which both sets of data can be reconciled. Kim and Lee []27[ analyse these two problems and 
present a solution based on maximally mixed( in vacuum) three generations of neutrinos. 
This later analysis is however a rather fine tuned solution since the parameter space allowed 



is rather tiny. 

In this paper, we analyze the solar neutrino problem by considering the oscillations 
between the three neutrino flavours. The analysis is done with no particular model of neu- 
trino masses and mixings assumed. The analysis is similar in spirit to that of Joshipura and 
Krastev |]26[| . We carry their analysis further and not only map out the full parameter space, 
but also discuss the average survival probability and recoil electron spectrum. In addition 
we also discuss a non-standard solution where no resonance occurs but nevertheless there 
is a parameter space in which all the three experiments discussed earlier can be reconciled. 
We also do not make any assumption about the evolution being adiabatic and take into 
account non-adiabatic effects. These effects may be ignored, however, in parts of allowed 
parameter space. In the three generation case the neutrino oscillations are determined by 
two mass differences and three mixing angles neglecting the CP-violating phase. One of 
the mixing angles is irrelevant for solar neutrino problem p3| , p4| while being relevant to the 
atmospheric neutrino problem and one of the mass differences is constrained by the atmo- 
spheric neutrino deficit. Therefore the solar neutrino oscillations in the three flavour case 
are dependent on three parameters. Because of the additional parameter, a larger region of 
the parameter space is allowed by the solar neutrino data compared to the two generation 
scenario. In section 2, we present the theoretical frame work for our analysis of the solar 
neutrino problem and in section 3, we present the numerical results for the solar neutrino 
problem in conjunction with the atmospheric neutrino problem. The last section consists of 
a brief summary and discussion. 



II. THREE NEUTRINO OSCILLATIONS IN MATTER- A PERTURB ATIVE 

ANALYSIS 

In this section we discuss the mixing between three flavours of neutrinos and obtain 
the probability for a v e produced in the sun to be detected as a v e on earth. The three 
flavour eigenstates are related to the three mass eigenstates in vacuum through a unitary 
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transformation, 











= U v 











(11) 



where the superscript v on r.h.s. stands for vacuum. The 3x3 unitary matrix U v can be 
parametrized by three Euler angles (u,<f>,if}) and a phase. The form of the unitary matrix 
can therefore be written in general as, 

U v = U phase x U 23 (4j) x UM) x U 12 (u), 



where U, 



is the mixing matrix between ith and jth mass eigenstates with the mixing 



angle It has been shown that the expression for electron neutrino survival probability, 
integrated over the time of emission and of absorption, is independent of the phase and the 



third Euler angle ip f23|j24}| . They can be set to zero without loss of generality and we have 
the following form for U v 



I \ 



JJV 







(12) 



where = sin0 and = cos0 etc. The angles uj and cf) can take values between and 7r/2. 
Note that one of the flavours decouples if either uo or is zero and we have a two flavour 
scenario. As mentioned earlier the approach here is similar to that of Joshipura and Krastev 
p6| who, however, assume that the mixing angle between second and third generation, ip, is 
small and hence can be neglected. We wish to emphasise that this is not an assumption and 
infact ij} can be arbitrary and the result for survival probability of the electron neutrino is 
independent of this |2^J2^]. In fact the solution of the atmospheric neutrino deficit requires 
ip to be rather large. Together, solutions of the atmospheric neutrino deficit and the solar 
neutrino problem determine the mixing matrix U v completely apart from the CP-violating 
phase. 
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The masses of the vacuum eigenstates are taken to be fix, fi% and fi%. In the mass 
eigenbasis, the (mass) 2 matrix is diagonal, 



Mi 



( fi\ ^ 
fi\ 

V A) 
( 



ft\i + 




6 2 i 

$31 



(13) 



where $21 = fi 2 ~ lA an d <^3i = lA ~ lA- Without loss of generality, we can take $21 and #31 
to be greater than zero. Neutrino oscillation amplitudes are independent of the first term 
so we drop it from further calculation. In the flavour basis the (mass) 2 matrix has the form 



M 2 = U v MlU v] 



°31 

M 3 x + 5 2 xM 21 , 



(14) 



where 



Mi 



31 



M, 



21 



s 







c\ J 



2 2 

C 4> S U1 C <f> S U} C U} 



Q 2 „2 

S <j> S u J 



(15) 



As in the two flavour case, matter effects can be included by adding A(r), defined in ( 
to the e — e element of M 2 The matter corrected (mass) 2 matrix in the flavour basis is 



Mi 



5 31 M 31 + 5 21 M 21 + AM A , 



(16) 



where 
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'lOO^ 



V 



(17) 



/ 



M A = OOO 


To calculate the evolution of a neutrino in matter we have to find the matter corrected 
eigenstates by diagonalizing M^. For arbitrary values of <5 3 i and 8 2 i, it is cumbersome to 
find the eigenvalues and eigenvectors of algebraically. However, the eigenvalue problem 
can be solved using perturbation theory, if the mass differences have the following hierarchy 
Ssi ^> 52i- This assumption is plausible in light of the observed atmospheric muon neutrino 
deficit. Recently Kamiokande analyzed their atmospheric neutrino data, assuming that the 
deficit is caused by the oscillation of a into another flavour. Their analysis assumes 
mixing between only two flavours [y^ <-> v e or «-> v T ). For both cases their best fit yields 
a mass square difference of the order of 10~ 2 eV 2 pT| . In our analysis we take <5 3 i to be 
10~ 2 eV 2 . Thus we have £31 much larger than A max and hence the oscillations involving 
the third generation are not influenced very much by the matter effects. In order for the 
matter effects to be significant (as necessitated by the solar neutrino problem), the other 
mass difference in the problem, 821, should be such that the resonance condition is satisfied 
for some values of parameters. This means S21 ~ A max . Thus we work in an approximation 
where 8 2 i,A max < 5 31 . 

In this approximation, to the zeroth order, both the matter term and the term pro- 
portional to 821 can be neglected in eq. fll6|) . Then = £ 31 M 31 , whose eigenvalues and 
eigenvectors are 

/ \ 

Ccj, 
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/ \ 

S31; • (18) 

Treating AM a + S21 M21 as perturbation to the dominant term in and carrying out 
degenerate perturbation theory, we get the matter dependent eigenvalues and eigenvectors, 

/ \ 

c<f> m c UJm 



m 



l! 



m\; 



m\; 



y s 4> m c u) m J 
( \ 

C 4>m S ^m 



y s 4> m s uj m J 
( \ 



C4, 



(19) 



The above eigenvectors are the columns of the unitary matrix U m which relates the flavour 
eigenstates to matter dependent mass eigenstates v™ through the relation 



(20) 



The matter dependent mixing angles can be expressed in terms of the vacuum parameters 
and A as 







V™ 




= U m 


m 
u 2 






m 
U 3 



tan 2uj r , 



621 sin 2oj 



1 + — COS 
031 



sin <p m = sin 
The matter dependent eigenvalues mf are given by 



821 cos Iuj — A cos 2 ' 
; cos m = cos 



1 A ■ 2 
1 — — sin . 

031 



(21) 
(22) 
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m\ = A cos 2 4> cos 2 u m + 821 sin 2 (u — uj m ) 
m\ = A cos 2 sin 2 u m + 8 2 i cos 2 (a; — a; m ) 
m 3 — ^31 + A sin 2 — ^31 • 



(23) 



u m can undergo a resonance if the values of #21, <p an d 00 are such that the resonance condition 



A{r) cos 2 (ft = 5 2 i cos 2a; 



(24) 



is satisfied for some r. Note that this condition is very similar to the resonance condition in 
the two flavour case (eq. The new feature here, which occurs due to the mixing among 
the three neutrino flavours, is the presence of the second mixing angle <ft m the resonance 
condition. This dependence on <fr leads to a larger region of allowed parameter space in the 
three flavour oscillation scenario as will be shown in the next section. Since S21, A(r) and 
cos 2 are all positive, a resonance can occur only if cos 2u is also positive, or if uo < 7r/4. 
In the three flavour case, the electron neutrino survival probability is given by 



[Pee) = E TO 



jjm 



(25) 



i>1 1 vj 1 ) is the probability that the jth matter dependent eigenstate evolves into zth 
vacuum eigenstate. As in the two flavour case, if the adiabatic approximation holds, then 



Sij. 



(26) 



We introduce the jump probabilities 



(z^K m > fori^j 



(27) 



to take into account the non-adiabatic transitions, if the adiabatic condition doesn't hold. 

Because £31 ^> A max ,62i, the third eigenvalue, both in vacuum and in matter, is much 
larger than the other two eigenvalues. Non-adiabatic effects are significant only if the eigen- 
values of two states come close together [ 28 1 . Therefore the jump probabilities involving the 
third state, X13 and X23 are expected to be negligibly small. Thus we have the expression 
for electron neutrino survival probability to be 
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- 2 6 

•rm 



{Pee) = cos 2 cos 2 m (cos 2 uj cos 2 u m + sin 2 uj sin 2 uj^ + sin 2 sin 2 

— %i2 cos 2 cos 2 (f) m cos 2tu cos 2cu m . (28) 

For X12 we use the formula, 

exp[-^]-exp[-^] 
" " l-exp[-^] ' (29) 

^ L 2 sin 2 u) J 

where 7 is defined in equation (§) and 

F = l-tan 2 cu (30) 

for an exponentially varying solar density |^| . We use this form for the jump probability 
since it is valid both for large and small mixing angles. In the extreme non-adiabatic limit 
X12 — > cos 2 uj and when jF » 1, we have the usual Landau-Zener jump probability given 
by X12 exp[— as expected. Infact for much of the allowed parameter space, this form 
can be used without any appreciable change in the results obtained. 



III. RESULTS 

In this section we discuss the results of the numerical analysis first for the solar neutrino 
problem and using that we map out the region in the parameter space which contains the 
solution to the atmospheric neutrino problem. 



A. solar neutrinos 

We analyze the expression for (P ee ) in (|28|) and find the ranges of 821, uj and <ft allowed 
by the three solar neutrino experiments. Since £31 ^> A max , we see from the expression for 
<p m in ( p2|) that the angle <fi is almost unaffected by the matter effects. However, u m can be 
significantly different from u and can undergo resonance if the resonance condition in ( p4|) 
is satisfied. Since this resonance condition depends on 0, in addition to 5 2 i and 00, a larger 
region of parameter space satisfies the three constraints from the experiments. 
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To search for the regions allowed in the three parameter space 821, to and 0, we define 
the suppression factors observed by the three types of experiments 

yGa = ^ Ga ^ = 0,563 ±0.067, 

R-Ga;SSM 

Vci= D RCI =0.318 ±0.051, 

VKam = D Kam = 0-51 ± 0.07, (31) 
J^Kam;SSM 

where the first number refers to the average of the data given by two experiments- namely 
GALLEX and SAGE. The predicted SSM rates for various experiments were taken from 
Bahcall-Pinsonneault SSM calculations 0. The uncertainties in y^ are the sum of the ex- 
perimental uncertainty in the numerator and the theoretical uncertainty in the denominator, 
added in quadrature. 

The predictions for yi for the three flavour oscillation scenario are obtained by convoluting 
the SSM fluxes and the detector cross sections with (P ee ) from The expression we use 
is 

= Ek It: x dE$ K (E)a(E) < P ee > (E) 
T.KSt-dE^ K {E)a{E) 

where the sum over K refers to the neutrino fluxes from various sources contributing to the 
process. We also include the contributions from the CNO cycle apart from the dominant 
contributions from the p-p cycle. In the case of Kamioka, only the 8 B flux contributes and 
one must also take into account the neutral current contribution arising from the muon 
neutrinos interacting with the detector material. The parameter ranges are then calculated 
by putting vetos on y at la and 1.6a levels. The energy dependent fluxes were taken from 
Ref. || and the cross sections were taken from Ref. P^] . 

Figure 1 shows the allowed values of u and <\> with 821 varying between 10 -6 eV 2 and 
10 -4 eV 2 . Note that the allowed values of 821 are also determined by the same veto con- 
ditions. In the two generation case it is a standard practice to plot 821 against sin 2 (2a;) 
since that is the combination that enters the survival probability. In the three generation 
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case all possible circular functions of the mixing angles are possible. Hence we depart from 
the standard practice in this paper and plot the angles themselves. The points refer to the 
allowed values after the vetos corresponding to all three experiments are imposed. The dark 
squares show the values allowed by la uncertainties given in (|3T| ) whereas the hollow squares 
show the values allowed when the uncertainty is increased to 1.6a. Fig.2 shows the allowed 
regions in the 0-521 plane, with u varying between and 7r/2 but obeying the same set of 
vetos. In Figs. 1 and 2 if we restrict ourselves to the = lines (the y-axes) we get the 
known |2| two-flavour solutions for u and £21- The large extended regions of the parameter 
space brought in through the additional degree of freedom in the three-flavour scenerio 
are shown clearly in Figs.l and 2. For completeness we also plot in Fig. 3 the allowed 
range in the ui — 821 plane. Here again the three-flavour scenerio provides an enlargement of 
the allowed parameter space over that of the two-flavour solution (small regions around the 
isolated dark patch in the left and around the end of the dark arm on the right). 
The various regions of the allowed parameter space may be classified as follows: 

1. small ^21, small u, small 0, 

2. large <5 2 1, large u>, small 0, 

3. small 621, small u, large 0, 

4. large S 2 i, small cu, large 0, 

5. large 62 1, large u, large 0, 

where the small or large S21 means either 621 < 10~ 5 eV 2 or 821 > 10~ 5 eV 2 . The first 
two regions corresponding to small in the above classification belong to an approximate 
two generation situation since the angle is small. The one corresponding to small us is the 
usual non-adiabatic solution, whereas the one corresponding to large u is the usual adiabatic 
solution. The rest invoke the genuine three generation oscillation mechanism. In the two 
flavour scenario, the small angle solution (corresponding to to small as in case 1 above) 
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gives the best fit @. There the parameter space allowed at la level is very small because 
the resonance condition and the non-adiabatic jump factor fix 5 2 i an d oj almost uniquely. 
These values of parameters indicate that the neutrinos from the p — p cycle suffer very little 
suppression and those from 7 Be suffer almost complete suppression as will be illustrated 
soon in the analysis of the survival probability. 

In the three flavour scenario, the resonance condition (eq. |4]) and the survival probability 
(eq. |28[) are dependent on the second angle 4> a\so. The suppression of the p — p neutrinos 
depends on the value of <fi and if this suppression is significant, then the complete suppression 
for 7 Be neutrinos can be relaxed. This is one of the important differences between the three 
flavour and the two flavour oscillations. 

Figure 4 shows the energy dependence of (P ee ) for some representative values of lo, 4> , 
and #2i- The curve labelled (a) corresponds to = 2°. As there is very little mixing between 
the first and the third generation of neutrinos, this is infact an almost two generation case. 
In agreement with the two generation analysis, there is almost no supression of the p — p 
neutrinos and the 7 Be neutrinos are almost completely suppressed. The survival probability 
at high energies relevant to Kamioka is almost a linear function with an average around 
0.5 as one would expect. Also here the values of uj and 821 are small (they are almost 
equal to the values obtained in the two flavour case) and the non-adiabatic effects become 
important beyond 2 MeV. Keeping uj small if we increase (f) in the allowed region there 
is a perceptible reduction in the probability in the p — p energy range and an increase in 
the survival probability of the 7 Be neutrinos (curves (b) and (c)). When 5 21 is increased, 
however, there is a qualitative change in the survival probability profile. In this range both 
uj and are allowed to be large. Here also there is a qualitative change when uj is small or 
large. For large uj the survival probability is a smooth function resembling the adiabatic case 
of the two generation analysis (curves (d) and (f)) whereas for small uj it is almost a step 
function (curve (e)) which is like the classic adiabatic case discussed by Bethe in the two 
generation case [|19|| . One common feature of the large 821 case is that the p — p neutrinos 
undergo substantial suppression varying between 0.6 -0.5. The resonance also occurs at a 
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much higher energy than in the small 821 case. Curve (f) has u, and 821 all large and 
in some sense it can be called 'most representative' of the three flavour oscillation scenario 
because both the mixing angles in this case are large. In all the above cases, except (e), 
the average survival probability above 7 MeV is in the neighbourhood of 0.4 which is what 
is required by the Kamioka data and there is no dramatic change from one to the other. 
This is not so at low energies where the curves differ dramatically. In this sense Kamioka 
experiment cannot distinguish between different theoretical scenarios of masses and mixings. 

One way of experimentally measuring the energy dependence of (P ee ) is to look at the 
recoil electron spectrum in those detectors that use u e — e scattering. In Fig. 5 we have 
shown the recoil electron spectrum for the six cases plotted in Fig. 4. Except case (f), they 
cannot be distinguished beyond 10 MeV, whereas there are substantial differences at low 
energies. While this energy range is not completely accessible in Kamioka, it is interesting 
to note that it may be possible to see this difference in the experimental recoil electron 
spectrum in the SNO JJIJ and Borexino |JT] detectors. Note that in computing the recoil 



electron spectrum, we have used the spectrum of 8 B neutrinos as input. This is because 
the threshold in experiments which can measure the recoil electron spectrum (like SNO 
and Kamioka) is more than a few MeVs where only this flux matters. The only exception is 
Borexino where the threshold is much lower and there are other contributions below 1.5MeV. 
In particular the 7 Be neutrino source, which is a line spectrum at 0.862 MeV, will show up as 
a sharp bump in the recoil spectrum where the height of the bump depends on the survival 
probability. A complete absense of the bump would point to the set parameters as in case 
(a) of Fig.4. 

Finally we consider a non-standard mixing which leads a substantial region in the pa- 
rameter space. We consider a situation where the electron neutrino is coupled more strongly 
to the heavier mass eigenstate v<i- Obviously this implies that the mixing angle between the 
first two mass eigenstates u, is greater than 7r/4. In the standard analysis the mixing has 
to be less than 7r/4 so that the resonance condition is satisfied as can be seen from eq. ([I4|) . 
This is true in the two as well as in the three generation case since the LHS of the resonance 
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condition is positive for arbitrary whereas the sign in RHS depends on the magnitude of u. 
There are no strong theoretical reasons not to consider this situation. Because u m at core is 
very close to ir/2, oj rn is then constrained to be ui < u) rn < n/2. Infact for 8 2 i < 10~ 7 eV 2 , u m 
is approximately n/2. Since hardly varies with density the effective survival probability 
may be approximately written as, 



where we have retained the jump probability xu- While it may appear some what unusual 
to keep the jump probability when there is no resonance, a plot of the eigenvalues clearly 
shows that the difference between the first two eigenvalues is not very different from that 
of the standard case close to vacuum and one cannot completely discard the existence of 
non-adiabatic jumps between mass eigenstates. However, most of the derivations of the 
jump probability assume the existence of resonance and the profile of the density variation 
close to resonance. Since we do not have a handle on this, we assume that the jump 
probability is simply given by rr 12 = exp(— C/E) and treat C as a free parameter of the 
theory. The survival probability is then energy dependent as would be required by the 
solution to the solar neutrino puzzle. The resulting parameter space is shown in Fig. 6 for 
u, <p. The parameter C varies from 0.4 to 6.3 in the allowed region. If we assume any one of 
the expressions for the jump probability discussed earlier, then we will have to discard small 
values of C (C < 4) since then the jump probability becomes very large and unacceptable. 
However in the allowed region, the points corresponding to small C are very few and there 
is no substantial change from the plot shown in Fig. 6. We also show some typical variation 
of the survival probability (P ee ) for some typical values of u, <fi and C in Fig. 7. The curve 
(a) corresponds to small C where non-adiabatic effects are important while the curve (b) 
corresponds to large C which is an almost adiabatic case. We wish to stress that this is an 
adhoc solution but we have analysed this situation because there are no strong theoretical 
reasons to ignore this possibility. 

To conclude this section, we note that the solution to the solar neutrino puzzle fixes 




(33) 
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the parameter space defined by u, <fi and 821 ■ While we have actually chosen the fourth 
parameter <5 31 to be 10~ 2 eV 2 we might as well have set the limit <5 31 > 10~ 3 eV 2 without 
affecting our results. One therefore requires more inputs to fix the range of 631 and the angle 
if) (mixing angle between second and the third generation neutrinos) which is arbitrary as 
far as the solar neutrino puzzle is concerned. The new input is provided by the analysis of 
the atmospheric neutrino problem which we consider next. 



B. Atmospheric Neutrinos 

In order to fix the mixing matrix completely we still need to fix the range of if), which 
is the mixing angle between the second and third generation neutrinos, as this is arbitrary 
in the solar neutrino analysis. To have a consistent solution for both solar neutrino and 
the atmospheric neutrino problems, we need to show that there exists a range of ip in the 
allowed range of parameters occuring in the solar neutrino problem. To ensure this we first 
define the ratio 

R = 4^ (34) 

which measures the ratio of the observed muon neutrino flux to the electron neutrino flux 
to that expected from Monte-Carlo calculations of neutrino production in the atmosphere. 
The most recent measurement of this ratio by the Kamiokande collaboration ]ll| yields 



R = 0.57^o;q7 ± 0.07 in the multi-GeV range. The depletion is further confirmed by the 
observation of the zenith-angle dependence. The result for sub-GeV range atmospheric 
neutrinos is R = 0.601^05 which is consistent with the multi-GeV range data. We may 
therefore assume that the suppression is approximately energy independent from sub-GeV to 
multi-GeV range of energies. Assuming that this depletion is due to the vacuum oscillations 
amongst the neutrino flavours, this ratio may be written as, 

P + rP 
P + -P 

1 ee T -i fie 
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where r = 4>v e /4> v = 0.45 is the ratio of the flux of electron neutrinos to that of muon neu- 
trinos at the production point. Note that r is simply the inverse of the flux ratio expected 
on the basis of the Monte-Carlo calculations (see eg. (|34|) ) . We now assume that the survival 
probabilities (P ee , P w ) and the oscillation probability (P^ e ) are given by the full three gen- 
eration mixing matrix defined by the angles (f), uj,ip and the two mass squared differences 
821, 831. The vacuum oscillation probability between two flavours is then given by, 

Pij =UfiUfi + Ul U% + U 2 3 U% + 2U a U i2 U jl U j2 cos(2.53-Ji) (36) 
+2U ii U a U jl U jZ cos(2.53-Ji) + 2U i3 U i2 U j3 U j2 cos(2.53-^), 

where i and j are the flavour indices, E is the energy given in units of MeV, 8ij is the mass 
differences in eV 2 and d is the distance of traversal given in meters. These probabilities 
explicitly depend on the distance d travelled by the neutrinos from the point of production 
to the point of detection and is approximately about 13,000 kms for the upward moving 
neutrinos. This distance is much less than the oscillation length between the first two 
generations( since 8 2 \ is small). Therefore the cosine factor involving 8 2 \ can be safely set 
equal to unity. As mentioned before, Kamioka has also observed that the level of suppression 
for the atmospheric muon neutrinos is approximately the same both for sub-GeV and the 
multi-GeV neutrinos. This can be ensured if the energy dependent factors involving 831 and 
8 32 are such that the cosine functions can be replaced by the corresponding averages. This 
is possible if and only if many oscillation lengths are contained in the distance travelled 
by neutrinos to the detector. This then sets the limits on the mass squared difference 
8 3 i > 10~ 3 eV 2 . The large 8 3 i regions (8 3 i > 10 _1 e^ 2 ) are excluded at 90 percent C.L by the 
analysis of the multi-GeV neutrino data . While we have used the central value 10~ 2 eV 2 
in our solar neutrino analysis, the results for both solar and atmospheric neutrinos remain 
unchanged if the value is further increased and marginal changes occur for values close to 
lO -3 ^ 2 because of the approximations we made in the solar neutrino analysis. 

Therefore the only range to be fixed is for the mixing angle ip. This we do by requiring 
the theoretical value of R calculated from eq.(^) is within la and 1.6a of the experimental 
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value. The resulting range for ip is shown in Fig. 8, whereas usual the full squares show the 
la veto and the open squares show 1.6a veto. 

A few comments are in order here: As in the two generation analysis of the atmospheric 
neutrino problem, we find that the preferred values of ip is large and around 7r/4. This can be 
checked easily by looking at the conversion probability P^ e in the allowed range of parameters 
for the atmospheric neutrino problem. It turns out that this conversion probability is always 
less than twenty percent. Thus the solution to the atmospheric neutrino problem is mainly 
driven by the — v T oscillations whereas the solution to the solar neutrino problem is mainly 
driven by u e — oscillations at least for small values of 0. However there are large domains 
of the parameter space where one requires the full three generation analysis presented here, 
to have a consistent solution to both the problems. 

IV. SUMMARY AND DISCUSSION 

We have examined in detail the possible solutions to the solar neutrino and atmospheric 
neutrino puzzles in the realistic three generation framework. There are in general three 
mixing angles, one phase from the mixing matrix and two mass squared differences which 
define the full parameter space. In the case of solar neutrinos the survival probability for 
the electron neutrino, even after taking into account the matter effects, is independent of 
the phase and one of the mixing angles. We also fix one of the mass squared differences 
by appealing to the atmospheric neutrino problem. Thus our parameter space in the so- 
lar neutrino analysis consists of two angles and one mass squared difference. In our case 
these are chosen to be u which gives the mixing between first and second generations , 
which is the mixing between first and third generations and #21 which is the mass squared 
difference between the first two generations. The mass difference is fixed to be around 
10~ 2 eV 2 to explain the atmospheric neutrino problem. We have mapped out the parameter 
space(0, u, 621) by invoking the vetos arising from the data given by the three solar neu- 
trino experiments. Next we have used these allowed ranges of parameters from the solar 
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neutrino analysis as input in the atmospheric neutrino analysis to fix the angle ip and find 
that there exists a substantial range in this parameter which allows a solution to the atmo- 
spheric neutrino puzzle. The numerical calculations necessarily depend on the bin size for 
the parameters. We have ensured that the bin size we have chosen is such that a further 
reduction will not change the overall profile of the allowed region. However it is conceivable 
that the rough edges that one still sees in parts of the allowed region will be smoothed out 
by a further reduction of the bin size. 

In conclusion, we have shown that there exists a consistent solution to the solar and 
atmospheric neutrino deficit puzzles within the framework of standard MSW mechanism 
based on the set of all available measurements of the solar neutrino fluxes. The full anal- 
ysis involves five parameters which we have mapped out by accommodating the solar and 
atmospheric neutrino fluxes seen by the present set of experiments. While the allowed re- 
gion in the parameter space is still large, these can be constrained further by measuring the 
distributions of recoil electron energies in solar neutrino detectors that use v — e scattering. 
Although the threshold energy at the Kamioka detector is rather too high for this purpose, 
the SNO and Borexino detectors may be effective in narrowing the parameter space. Fi- 
nally we would like to remark that the analysis of solar and atmospheric neutrino problems 
presented here is exploratory in nature. This is so since with time the errors are bound to 
change which inturn will affect the vetos imposed by us at la and 1.6a levels. Nevertheless 
we believe there is already sufficient indication that a robust solution of both problems is 
possible within the framework provided by the mechanism of neutrino oscillations with three 
generations. 
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FIGURES 

FIG. 1. Allowed regions in (j> — u plane (with 10~ 6 eV 2 < 821 < 10~ 4 eV 2 ) at la (dark squares) 
and at 1.6<r (hollow squares). 

FIG. 2. Allowed regions in — log(<52i/eV 2 ) plane (with < ui < vr/2) at la (dark squares) 
and at 1.6<r (hollow squares). 

FIG. 3. Allowed regions in w- log((52i/eV 2 ) plane (with < < ir/2) at la (dark squares) 
and at 1.6a (hollow squares). 

FIG. 4. Survival probability (P ee } vs E v for typical values of 0, uj and 821 in the al- 
lowed region. The parameters chosen are: (a) 621 = 4.0 x 10~ 6 ,w = 2.5°, <p = 2.0°; (b) 
5 21 = 5.0 x 10~ 6 ,u; = 2.0°, (f> = 16.5°; (c) 5 2 i = 7.0 x 10~ 6 ,lj = 1.75°, <j> = 37.5°; (d) 
621 = 2.5 x 10~ 5 ,lj = 35.0°, (p = 3.0°; (e) S 2 i = 7.0 x lO' 5 , uj = 2.0°,^ = 30.0°; (f) 
621 = 1.0 x 10" 4 ,cj = 24.5°, = 24.0°; 5 2 i is given in terms of eV 2 . 

FIG. 5. Recoil electron spectrum for different representative points of the allowed parameter 
region. The parameters for the differenct curves labelled (a)-(f ) are the same as in Fig. 4. The inset 
shows a comparison of all zix cases with the SSM spectrum(dashed line). 

FIG. 6. Allowed regions in <p — uj plane (with0.4 < C < 6.4 (dark squares) and at 1.6<t (hollow 
squares) for the non-standard solutions. 

FIG. 7. Typical survival probability profile in the non-standard case. The curve la- 
belled (a) corresponds to C = 0.4, u = 55°, 4> = 2° and the curve labelled (b) corresponds to 
C7 = 6,cj = 89°,(/) = 38 . 

FIG. 8. The allowed range of values for the mixing angle i/j in the i/j — <p plane when the (f> and 
uj are restricted to the range allowed by the solar neutrino problem(see Fig.l). 



25 



